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ABSTRACT 

Context. The formation of massive stars is a highly complex process in which it is unclear whether the star-forming gas is in global 
gravitational collapse or an equilibrium state supported by turbulence and/or magnetic fields. 

Aims. By studying one of the most massive and dense star-forming regions in the Galaxy at a distance of less than 3 kpc, i.e. 
the filament containing the well-known sources DR21 and DR21(OH), we attempt to obtain observational evidence to help us to 
discriminate between these two views. 

Methods. We use molecular line data from our '^^CO 1— >0, CS 2— >1, and NiH^ 1— >0 survey of the Cygnus X region obtained with 
the FCRAO and CO, CS, HCO+, N2H+, and HjCO data obtained with the IRAM 30m telescope. 

Results. We observe a complex velocity field and velocity dispersion in the DR21 filament in which regions of the highest column- 
density, i.e., dense cores, have a lower velocity dispersion than the surrounding gas and velocity gradients that are not (only) due to 
rotation. Infall signatures in optically thick line profiles of HCO^ and '-CO are observed along and across the whole DR21 filament. 
By modelling the observed spectra, we obtain a typical infall speed of ~0.6 km s"' and mass accretion rates of the order of a few 
10"^^ Mq yr"' for the two main clumps constituting the filament. These massive clumps (4900 and 3300 Mq at densities of around 
10^ cm"^^ within 1 pc diameter) are both gravitationally contracting. The more massive of the clumps, DR21(0H), is connected to a 
sub-filament, apparently 'falling' onto the clump. This filament runs parallel to the magnetic field. 

Conclusions. All observed kinematic features in the DR21 filament (velocity field, velocity dispersion, and infall), its filamentary 
morphology, and the existence of (a) sub-filament(s) can be explained if the DR21 filament was formed by the convergence of flows 
on large scales and is now in a state of global gravitational collapse. Whether this convergence of flows originated from self-gravity on 
larger scales or from other processes cannot be determined by the present study. The observed velocity field and velocity dispersion are 
consistent with results from (magneto)-hydrodynamic simulations where the cores lie at the stagnation points of convergent turbulent 
flows. 
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1. Introduction 

Our goal in this series of papers is to investigate the generic link 
between the large-scale (several 10 pc) molecular-cloud spatial 
structure (Schneider et al. 2010), medium-scale (<10 pc) dy- 
namic fragment properties (this paper), and the occurrence of 
high-mass star formation on small (0. 1-0.5 pc) scales (Bontemps 
et al. 12010 Csengeri et al. ,2010) . To achieve this, we make use 
of our multiwavelength study in Cygnus X (see below). This re- 
gion has akeady been shown to represent an excellent laboratory 
for studies of high-mass star formation. 

Molecular clouds (MCs) are the birthplaces of low- and high- 
mass stars. Understanding their formation and evolution is es- 
sential to understanding in general the star formation (SF) pro- 
cess. There are currently two different scenarios for the evolution 
of MCs and stars, a dynamic and a quasi-static view. 

In the dynamic context, MCs are transient, rapidly evolving 
entities that are not in equilibrium, and the spatial and velocity 
structure of the cloud is determined by compressible supersonic 
turbulence (see e.g., Mac Low & Klessen l2004l for an overview). 
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The driving sources of turbulence can be diverse and their rela- 
tive importance remains a subject of debate. Large-scale mecha- 
nisms such as supernovae explosions, should be the most impor- 
tantXS, based on theoretical arguments extensively discussed by 
MacLow & Klessen ([2004] |. 

Turbulence may also occur during the formation process of 
molecular clouds (Vazquez-Semandeni et al. 2002, Klessen & 
Hennebelle 12010b within large-scale colliding flows of mostly 
atomic gas in the galactic disk, generated by dynamic compres- 
sion in the interstellar medium or other large-scale instabilities 
(see, e.g., Hennebelle & Audit [2007, Vazquez-Semadeni et al. 
2008, Hennebelle et aL l2008l Banerjee et aL|2009l). High veloc- 
ity compressive flows form dense structures at stagnation points 
that may collapse to form stars/clusters (Ballesteros-Paredes et 
al. .200 3 , Vazquez-Semadeni et al. 2007, 2009 Heitsch et al. 
120081) . In this case, the velocity field and velocity dispersion 
of molecular clumps may contain the signature of the external 
convergent motion and the compressive, gravitational contrac- 
tion (this was already proposed by Goldreich & Kwan fT974] l. 
The star-forming core is then the dense post-shock region with 
a more quiescent velocity dispersion than the turbulent flow 
(Klessen et al. 12005b . This scenario is valid for both, low- 
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and high-mass star formation. However, there is no direct ob- 
servational confirmation of the existence of convergent flows. 
Studying the dynamics of HI and molecular gas may be a possi- 
bility (Brunt 2003 ), though it is not clear what observational sig- 
natures are to be expected. An indirect argument for molecular 
cloud formation out of colliding HI flows was provided by Audit 
& Hennebelle (2009 ), who showed that their models closely re- 
produce the observed clump mass spectra and Larson-relations. 
It has also been noted that the lifetime of the cloud in the gravo- 
turbulent framework is short (one dynamical crossing time, i.e. 
~10^ yr for giant molecular clouds) and the star formation pro- 
cess is rapid. 

In the quasi-static view (see McKee & Ostriker l2007l and ref- 
erences therein), MCs are formed by large-scale self-gravitating 
instabilities such as spiral density waves. (High-mass) star for- 
mation is approached in the 'turbulent core' model (McKee & 
Tan 2003) in which the star-forming cores are supersonically 
turbulent. Turbulent magnetic and thermal pressure supports the 
clump against self-gravity. Because of energy injection from 
newly formed stars, the clumps and most of the cores may main- 
tain their equilibrium, unless they quasistatically evolve to a 
gravitationally unstable state to finally form stars (if the mag- 
netic flux diffuses out of the clump by ambipolar diffusion). In 
the case of low-mass stars, the collapse of a rotating cloud of 
gas and dust leads to the formation of an accretion disk through 
which matter is channeled onto a central protostar For stars with 
masses higher than about 8 M© this mechanism of star forma- 
tion is less clear due to the strong radiation field that pushes 
against infalling material and may halt accretion. However, the- 
oretical work (Yorke & Sonnhalter 2002; Krumholz, McKee, 
Klein 12005; Peters et al. 2010) has shown that outflows lead to 
anisotropy in the stellar radiation field, which reduce the radia- 
tion effects experienced by gas in the infalling envelope. Thus, 
massive stars may therefore be able to form by a mechanism 
similar to that by which low mass stars form. Peters et al. ( 12010b 
even show that clustered SF is a natural outcome of massive SF 
even in the presence of radiative feedback. However, clustered 
SF is also a natural result of the competitive accretion scenario 
(Bonnell & Bate 2006 ), in which the fragmentation of a turbu- 
lent cloud produces stars with masses of the order of the Jeans 
mass within a common gravitational potential. These stars, lo- 
cated close to the centre of the potential, accrete at much higher 
rates than isolated stars and become massive. 

Both scenarios have in common that massive (>10Mq) stars 
form only from the coldest, and densest, cores (size scale <0. 1 
pc) of molecular clouds and that high infall/accretion rates are 
required to overcome feedback processes (ionizing radiation, 
jets, outflows). In the dynamic framework, stars form as the dy- 
namical evolution of the MC progresses with gas continuously 
being accreted. Individual clumps (~0. 1-0.5 pc) fragment into 
hundreds of protostars, possibly competing for mass in the cen- 
tral regions of the cluster (Bonnell & Bate l2006] l. Support from 
magnetic fields (Hennebelle & Teyssier 2008) and/or radiation 
(Krumholz 120061 Bate .2009) could drastically Hmit the level of 
fragmentation and channel the global infall to fewer, more mas- 
sive protostars (Bontemps et al.'2010'). 

From the observational point of view, indications of global 
collapse have been detected using molecular lines. Spectral pro- 
files of high density tracers, usually combining at least an op- 
tically thin and an optically thick line, are good probes of in- 
falling gas (e.g. Myers et al. 11996b . Observations in low-mass 
star-forming regions (e.g. Lee et al. p003 ) contain a mixture of 
infall (blue-shifted emission and/or redshifted self-absorption in 
the optically thick line) and outflow asymmetry. This sort of line 



profile was also found in the high-mass star-forming region W43 
(Motte et al. '2003). Peretto et al. ( 12006b proposed that a pecu- 
liar velocity discontinuity could be the result of some large-scale 
motion, originating from self-gravity (Peretto et al. 2007 ). The 
theory of turbulent core formation, on the other hand, predicts 
the existence of massive pre-stellar cores that have not been de- 
tected so far, not even in sensitive and extensive dust continuum 
studies (e.g. Motte et al. 2007 ). 

In this paper, we shortly introduce the Cygnus X region and 
the DR21 filament in Sect. |2]and describe the molecular line 
oberservations in Sect. [3] and show maps and spectra in Sect. 4] 
Section 5]presents an analysis of the kinematic structure and the 
physica properties of the DR21 filament and in Sect. [6] we use 
our findings to test the conditions of the different high-mass star 
formation and turbulence models. Section|7]summarizes the pa- 
per. 



Table 1. Observing parameters of the molecular line data ob- 
tained with the FCRAO and IRAM 30 m telescope: Column 
one and two indicate the molecular transition and frequency, fol- 
lowed by the half power beam width (HPBW) in arcsec, rj^^ is 
the main beam efficiency, Tsys the system temperature, Avres de- 
notes the velocity resolution, and AT^^s the average rms noise 
temperature per channel on a scale. 
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2. Cygnus X 

The Cygnus X region 

Cygnus X is one of the richest star-formation regions in the 
Galaxy, covering an area of about 30 square degrees in the 
Galactic plane around Galactic longitude 80° (Reipurth & 
Schneider 12008b . It contains a prominent OB-association, Cyg 
OB2, with ~100 O-stars and a total stellar mass of up to 10^ 
Mo (Knodlseder 2000) . From lai-ge-scale '^CO 2^1 (KOSMAQ 
Schneider et al. l2006l) and '^CO 1^0 (FCRACQ Schneider et 
al. 120071 Simon et al., in prep.) surveys, we obtained a mass of 
a few 10^ Mq for the whole molecular cloud complex that is di- 
vided into two parts - Cygnus North and South. These studies 
also showed that the majority of molecular clouds is located at a 
common distance of about 1.7 kpc, i.e. the distance of Cyg OB2. 

Its proximity makes Cygnus X one of the rare laboratories 
where different phases of massive star formation can be studied 
in detail. More than 800 distinct HII regions, a number of 

' Cologne Observatory for Submm-Astronomy 
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Fig. 1. Zeroth moment map of CS 2^1 emission (red contours, mapping area indicated by a continuous blue line) in Cygnus X 
overlaid on a moment map of '^'CO 1— >0 emission (mapping area indicated by a dashed line). Both maps were obtained with the 
FCRAO and have an angular resolution of ~50" . Emission was integrated in the velocity range -10 to 20 km s"' and the grey scaling 
ranges from 1 to 15 K km s Red contour levels for CS are at 0.5,1,2,4 K km s The DR21 filament is indicated, as well as some 
prominent objects in Cygnus X (S106, AFGL2591, W75N, the most massive stars of the OB2 cluster). 



Wolf-Rayet stars, several OB associations, and at least one 
star of spectral type 04If are known in Cygnus X, reflecting 
its record of high-mass star formation in the past. Ongoing 
massive star formation is revealed by wide-field (3 deg^) 1.2 
mm continuum imaging of the densest regions in the Cygnus 
X molecular clouds (Motte et al. 120071) . More than 100 proto- 
stellar dense-cores were detected of which 40 are likely to be 
the precursors of massive OB stars. High-angular resolution 
observations with the Plateau de Bure Interferometer (Bontemps 
et al. 1201 Oi l have indeed confirmed that the most massive dense 
cores of this sample are the sites of massive star formation, 
mainly in the form of clusters. 

The DR21 filament 

From our 1.2 mm continuum imaging and the '■'CO sur- 
veys, we concluded that the molecular ridge containing DR21 



and DR21(OH) (Dickel, Dickel & Wilson [T9781 Wilson & 
Mauersberger 1990, Jakob et al. 2007) is the most active, dense 
(average density~10'* cm"-'), and massive (34 000 Mg) cloud in 
Cygnus X. 

DR21 (Downes & Rinehart 1 19661 1 itself is a group of sev- 
eral compact HII regions (Harris et al. 11973b with an associated, 
very energetic outflow (e.g. Garden et al. I1991al Russell et al. 
1992), located in the southern part of the ridge. Three arcmin- 
utes further north lies DR21(OH), famous for its maser emis- 
sion (H2O, Genzel & Dow nes \M 2\ OH, Norris et al. fT982l 
CH3OH, B atria & Menten 19881. Even further norfli lies the 
massive star-forming region W75N (see Shepherd et al. l2004l for 
a review). Very recently, the DR21 filament was again the target 
for a wealth of studies in various wavelengths. Infrared images 
at 3.6, 4.5, and 8 fim from the Spitzer satellite (Marston et al. 
120041 Hora et al. 120091 ) show the complexity of the region with 
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Ri(aOOO) BA(3000) 

Fig. 2. Zeroth moment maps of CS 2^1 (left, contour levels are 0.5, 1, 2, 4 K km s ') and N2H"^ 1^0 (right, contour levels are 
0.5, 1.5, 3 K km s ') emission in Cygnus X North in the velocity range -10 to 20 km s ' observed with the FCRAO. Thermal 
HII regions (DR17-23) are indicated by red stars, mm-continuum sources from Motte et al. (2007) by green triangles. The latter 
correspond very well with peaks of NaH^ emission. 



many IR-filaments perpendicular to the ridge that seem to be 
sites of star formation (Kumar et al.|2()07 ). Dust continuum stud- 
ies (Vallee & Fiege 2006, Motte et al. 2007 ) confirmed the detec- 
tion of 8 compact, dense cores in the DR21 filament (Chandler et 
al.[l993 ) and found additional ones. Some of them are high mass 
protostars and drive outflows detected in SiO (Motte et al. l2007l l. 
An H2 1-0 S(l) image at 2.121 /jm (Davis et al. l2007] l shows at 
least 50 individual outflows driven by embedded low-mass stars. 

We studied this particular filament in high-angular resolu- 
tion observations with the IRAlvj^ 30 m telescope in different 
molecular tracers (isotopomeric CO lines, CS, C^'^S, HCO^, 
H'^CO+, N2H+, H2CO) to investigate the distribution of the 
different phases of molecular gas (cold dense cores, warm en- 
velopes) and to uncover infall/outflow signatures. The lower an- 
gular resolution FCRAO data delineate the large-scale structure 
in which the DR21 filament is embedded. A more detailed anal- 
ysis of the FCRAO '^CO data set is presented in Schneider et 
al. (2007 2010) and the CS and N2H+ data will be discussed in 
a forthcoming paper. All molecular line data we obtained serve 
as a basis for a physical model of the DR2 1 filament, including 
detailed non-LTE line modelling. 

3. Observations 

3.1. FCRAO 

We used data from a molecular line survey ('^CO and C'^O 
J=1^0, CS J=2^1, and N2H+ J=1^0) of the entire Cygnus X 
region (Simon et al., in prep., Schneider et al. i2007i) taken with 



the FCRAO 14m telescope. The data were obtained between 
2003 December and 2006 January using the 32 pixel 'Second 
Quabbin Optical Imaging Array' (SEQUOIA) in an On-the-Fly 
(OTF) observing mode. In this paper, we employ '^^CO 1— >0, 
CS 2^1, and N2H^ 1— >0 data that cover an area of 35 square 
degrees in '^CO and 12 square degrees in CS/N2H^ on a 22". 5 
grid. The beamwidth of the FCRAO at 93 (1 10) GHz is 45"(48") 
and the dual channel correlator was configured to obtain a veloc- 
ity resolution of ~0.1 km s"'. The spectra have a mean 1 cr rms 
noise level of ~0.2 K on a T^ antenna temperature scale, i.e., 
not corrected for the main beam efficiency of ~0.48. Pointing 
and calibration were checked regularly at the start of the Cygnus 
observing interval and after transit (no observations were per- 
formed at elevations higher than 75°). Pointing sources were SiO 
masers of evolved stars, i.e.,^-Cyg, R-Leo and T-Cep, depend- 
ing on LST-time. The calibration was checked regularly on the 
position of peak emission in DR21 and found to be consistent to 
within 10%. 

3.2. IRAM 

All molecular line map^ obtained with the IRAM 30m tele- 
scope in the Sierra Nevada/Spain were performed in June 2007. 
Table [T] gives an overview of all lines and observational param- 
eters. In order to efficiently map an area of ~4'xl0', covering 
the whole DR21 filament, we observed in the 'on-the-fly' map- 
ping mode. The total area was divided into 7 sub-maps that were 
scanned horizontally. For strong lines (e.g. CO, HCO^), only 
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The IRAM data are available on request, contact nschneid@cea.fr. 
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Fig. 3. Channel maps of '^CO 1— »0 emission in Cygnus X North in the velocity range -3.2 to -0.9 km s ' observed with the FCRAO. 
The red triangle shows the position of DR21(OH). The three major subfilaments (Fl, F2, F3) linked with the large DR21 filament 
are indicated in the plot of velocity -2.6 km s"'. 



one or two coverages were required while weaker lines were ob- 
served more often or with two receivers at once to obtain a suf- 
ficient signal-to-noise ratio. The resulting regridded (5") maps 
show slight smearing effects and fluctuations in the intensity cal- 
ibration only in the lines of '^CO and H2CO. The center po- 
sition of all maps is within 3" of the MM2 -position given in 
Mangum, Wooten, and Mundy (1991): RA(2000)=20''39'"00' 
and DEC(2000)=42°22'44". As OFF-position, we used a po- 
sition with ofi^set 0",1800" from DR21(OH) where '^CO 2^1 
emission of around 1 K is observed at velocities around -42 km 
s~\ separated well from the bulk emission of the DR21 filament. 
All data were reduced using GILDASj^ 

The spectral lines were recorded simultaneously in uniform, 
good weather conditions (average atmospheric opacity of 0.2 at 
230 GHz) prevailing over several days. We used the IRAM facil- 
ity receivers and autocorrelator, which were adapted in terms of 
spectral resolution and bandwidth such that all lines have a sim- 
ilar velocity resolution between 0.03 and 0.05 km s"'. Pointing 
and focus were checked every 2 hours. The pointing accuracy 
was found to be better than 4" and the receivers were aligned to 
within 2". 



4. Results 

4. 1 . The DR21 filament as part of a complex network 

Figure [T] shows the zeroth moment map of line integrated (-10 
to 20 km s ') CS 2— >1 emission obtained with the FCRAO over- 
laid on a grey scale map of '^CO 1^0 emission. The moment 

^ http:iram.fr/IRAMFR/GILDAS 



maps were constructed following the masked moment procedure 
described in Adler et al. (119921) . For CS and N2H+, we used 
a temperature threshold of 0.1 K, for '^CO a threshold of 0.2 
K. For the CS mapping, we focussed on the high column den- 
sity regions seen in '^CO. We adopt the nomenclature 'Cygnus- 
X North, South, and Northwest', established in Schneider et al. 
(120061 ) for the different cloud regions. 

The CS emission is more compact than that of '^CO and 
clearly delineates the high density regions. Interestingly, regions 
of peak emission in '^'CO do not always show up in CS. This 
indicates that, though the column density is high, the average 
density in these regions may stay below the critical density (n^r) 
of the CS line (a few 10^ cm"^^). Abundance variations caused 
by chemical effects probably also play a role. While in Cygnus 
North, and partly in its north-west region, CS clumps form a 
connected network of filamenatry structures, the emission distri- 
bution in Cygnus South is more dispersed with more isolated CS 
clumps. 

In the following, we concentrate on a discussion of the fil- 
aments in Cygnus-North. For that, we display in Fig. |2] maps 
of CS 2— >1 and N2H^ 1— >0 emission where the 'organization' 
of clumps along filaments is clearly visible. Even the NaH^ 
map - which traces only dense (n„-=0.5-70xl0*' cm"^) and cold 
gas (e.g. Tafalla et al. 2002) - shows the same arrangement of 
clumps. Most prominent in all maps is the filament containing 
the HII region DR21 (see below). The moment maps (Fig. |2|) 
show that all evolved star formation sites in Cygnus North (indi- 
cated by the HII regions DR17 to DR23) as well as future star- 
forming sites (indicated by the mm-continuum sources tracing 
massive dense cores, Motte et al. 120071 ) lie on filamentary struc- 
tures. 
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4.2. The DR21 filament 

From the large-scale, medium-angular-resolution (~1') molec- 
ular line maps (Sect. |4.1| i and previous studies (Sect. [TJ, the 
DR21 filament clearly appears to be the most prominent star- 
forming region in Cygnus X. The DR21 region in the south con- 
tains a very powerful, evolved outflow, while further north, the 
DR21(OH) region contains maser sources, indicative of recent 
star-formation activities. The mm-contiuum survey of Motte et 
al. (2007 ) revealed more than 20 protostellar objects, i.e. sites of 
star formation. The filament is shown in more detail in chan- 
nel maps of '^CO 1— »0 emission (Fig. [sjl, obtained with the 
FCRAO. It reveals the characteristic N-S oriented ridge known 
from other molecular line and continuum surveys (Dickel et al. 
[Mil Wilso n & Mauersb erger 1990, Vallee & Fiege 2006 Motte 
et al. 120071 Jakob et al. I2007I I. What is special in these chan- 
nel maps is the appearance of several E-W oriented subfilaments 
(indicated as F1-F3 at v=-2.6 km s"') that are connected to the 
main DR21 filament and extend over several velocity channels. 
The most prominent one (F3) has a length of around 10' (5 pc) 
and is also visible in the CS map (Fig|2|l and in mm-continuum 
(Motte et al. 12007) . indicating that gas of high density is present. 
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Fig. 4. Positionally averaged (across the whole observing region 
for each line) spectra of all lines observed in the DR21 filament 
using the IRAM 30m telesope. The HCO^ line is multiplied by 
a factor 8.5, the C^'^S and U^^CO+ lines by 10. 



4.2.1. Average spectra 

To study in more detail the physical properties (mass, density, 
temperature) and kinematics (infall and outflow signatures, ve- 
locity gradient) of the filament, we mapped the central part of 
the filament containing DR21 and DR21(OH) in various molec- 
ular line tracers at an angular resolution of between 12" and 30" 
using the IRAM 30m telescope. 

To distinguish the different emission features of the filament, 
we produced positionally averaged spectra of five of the ob- 
served lines (Fig. |4|i. The averaging was performed across each 
map extent, i.e. covering the whole filament. The CO lines ex- 
hibit two velocity components, one centered on -3 km s ' and 
one on +9 km s They correspond to the molecular clouds as- 
sociated with DR21 and W75N, respectively (see also Dickel et 
al. 119781 1. Though the W75N region itself is not included in the 
mapping region, molecular gas related to this star-forming re- 
gion northwest of the DR21 filament (see Fig.|2]i is still present 
across the filament. However, the gas is more diffuse and of 



lower density since a line at +9 km s"' (Fig.|4| can only be seen 
in the low-density tracers ('^CO and '^CO). The HCO+ line re- 
veals a narrow absorption feature at -1-9 km s"' (seen also by 
Nyman 1983), while the optically thin lines show no significant 
emission for the W75N component. 

The bulk emission of the DR21 filament is found at -3 km 
s ' (single Gaussian fines in C^'^S, H'^^CO+, and '^CO). The op- 
tically thick lines of '^CO and HCO^ show prominent, extended 
blue and red wing emission. The red wing, however, is blended 
with the W75N component, seen in emission in '^CO and ab- 
sorption in HCO^. The wings are caused by the powerful outflow 
ofDR21 (Gard en et al. 1991a) but also other outflow sources dis- 
The self-absorption feature at -3 km s"' is 



4.3 



cussed in Sect. 

seen in the optically thick lines of '^CO and HCO^, while opti- 
cally thin lines peak at the velocity of the absorption dip. This 
may be a signature of inflowing gas and is discussed in more 
detail in Sect. I 



4.2.2. The DR21 filament in velocity-integrated maps 

Based on the average spectrum (Fig. |4]) we produced line inte- 
grated maps of the bulk emission of the DR21 cloud (Fig. |5]). 
Though there are differences in the detailed emission distribu- 
tion, one can separate three main regions: a southern cloud with 
DR21, a middle emission peak with DR21(OH), and a northern 
extension with two secondary peaks. 

While DR21 is most prominent in the HCO^ and '^CO lines, 
DR21(OH) is more pronounced in all optically thin lines. A 
small shift in the position of peak emission of DR21(OH) can 
be discerned by comparing C^'*S/H'^CO+ and N2H+. For N2H+, 
we determined the optical depth using the hyperfine structure 
pattern of this line and found that t remains below 1 basically 
everywhere in the map. The latter traces cold, dense gas that is 
located more south of the warmer central part of the DR21(0H) 
region. Mauersberger et al. ( |1986t report a temperature above 
200 K in a small region of size <9" in DR21(OH), and Wilson 
& Mauersberger ( 119901 ) determine a rotation temperature of 34 
K from NH3 observations. Mangum et al. ( 11991119921 ) inter- 
preted their observations as a detection of a B-star within the 
DR21(OH) clump. The northern part of the filament is more 
prominent in NaH^ and H'^CO^ emission and not in C^'^S. 
Wilson & Mauersberger (11990b . and Wienen (120081) derive val- 
ues around 20 K for this region. 

Figure [6]is a position-velocity cut in '^CO and '^CO 2— >1 
and HCO^ 1— >0 along the filament in which all data points in RA 
were averaged. The greyscale image of '^CO emission shows 
that the bulk emission of the DR21 filament at -3 km s ' and the 
bulk emission of the W75N component at 9 km s"' are appar- 
ently connected. The gas linking the two components has a low 
column density since it is only visible in '^CO and not in '^CO 
(red contours). A lack of '^CO emission exists between A6=Q" 
and -150" at ~8 km s"' , indicating that the (cooler) W75N cloud 
lies in front of the DR21 filament and causes the '^CO self- 
absorption. Dickel et al. (1978) and Bieging et al. (119821) came 
to the same conclusion following their detection of self-absorbed 
line profiles in H2CO and OH toward DR21. However, the situa- 
tion is slightly more complicated because the W75N line compo- 
nent is blended in its southern part (<0") with emission from the 
'Great Cygnus rift', which is basically not seen in '^CO emis- 
sion (Schneider et al. 120071 ). The rift is an extended, low-density 
gas cloud at a distance of around 600 pc. Its emission can start 
as low as +4 km s"' and makes it thus difficult to disentangle the 
different emission features. It is, however, possible that there is a 
real physical connection/interaction between the DR21 filament 
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Fig. 5. Velocity-integrated spectral line maps of the DR21 filament. The HCO^ 1— >0 and '^CO 2— >1 lines are integrated over a 
velocity range from -11 to 5 km s"', the H2CO 3(1,2)— >2(1,1) line from -11 to 2 km s"', and all other lines from -7 to 1 km s"'. 
Contour levels follow the notation start/end/step with 15o-=3.8/87.3/7.5 K km s"' (HCO+), 6cr=l. 3/10.9/3.9 K km s"' (H'^CO+), 
30cr=34.1/249.4/23.9 Kkm s-i ('^CO), 6cr=5.9/14.7/2.9 K km s"' (C^'^S), 9cr=6.4/49. 1/10.7 K km s"' (H2CO), 3cr=0.27/l .9/0.27 
K km s ' (N2H+). Triangles indicate the positions of DR21(OH) and DR21. 
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Fig. 6. Position-velocity cut along constant RA, i.e. along the DR21 filament. Left: '^CO 2— >1 (grey scale and grey contours) and 
'^CO 2^1 (red contours). Contour levels are 2 to 34 K km s ' by 2 K km s ' for '^CO and 1.5 to 16.5 K km s"^ by 3 K km s"' for 
'^CO. The stars mark DR21(OH) and DR21. Right: The same plot for HCO+ with contour levels from 1 to 8 K km s^^ by 0.5 K 

km s"'. 



and the W75N cloud, as proposed by Dickel et al. (1978). They 
suggested that the 'W75N' cloud is differentially slowed down 
as it passes the edge of the DR21 cloud due to the observed N-S 
oriented velocity gradient from -nil km s"' to +9 km s"' (which 
is clearly seen in our CO data as well). 



Several outflows are identified out of which the DR21 out- 
flow at A^=-200" is the most pronounced. However, its red flow 
is not clearly defined because (i) it is blended with the 'Great 
Rift' emission and (ii) because it runs into an HII region (Garden 
et al. I1992I I. The blue flow, however, is clearly defined. Other 
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'flow-features' are identified at Ad-0" (DR210H(M)) and at 
~180", both look well-defined at negative velocities and dis- 
turbed for positive velocities. The position-velocity diagram for 
HCO^ shows that these outflows are much less obvious than that 
towards DR21. We come back to outflows in Sect. 14.31 and in 



particular the one seen at DR210H(M) in Sects. 4.2.3 and 5.2 



Another feature visible in Fig.|6]is a region of self-absorption 
in '^CO (emission regions in '^CO) at around -3 km s ' along 
the whole filament (indicated by a dashed white line). This cen- 
tral region of self-absorption is even more clearly identified in 
the HCO^ map. The strongest self-absorption is found at the 
dechnation off'sets of DR21 and DR210H(M). We discuss this 
finding in more detail in Sect. |5.2| 



channel maps of optically 



4.2.3. A subfilament at km/s ■ 
thin lines 



Figure|7]shows channel maps of H^^CO^ and C-^'^S emission be- 
tween + 1 and -6 km s The distribution of bulk emission of 
the cloud between -2 and -5 km s"' is similar in both tracers and 
follows the typical NS-ridge of the DR21 filament. Some emis- 
sion peaks are stronger in H'^CO^ than in C^'*S (and vice versa) 
due to the different density and temperature regimes. An exam- 
ple is seen at velocity -4 km s"' where C^"^S peaks prominently 



at the DR21(OH) clump while H'^CO^ emits more strongly in 
the northern part of the filament. 
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Fig. 8. Spectra of several line tracers at the position -25",25". 
The '^CO line was reduced by a factor 3.5 and for the N2H^ line, 
we only show the 101-012 component of the hyperfine structure 
(to avoid confusion). The different emission components are in- 
dicated with short dashed lines. 

At velocities larger than -2 km s"\ we identify an extended 
emission component between -1 and 0.5 km s"' that appears 



N. Schneider et al.: DR21 filament 



9 



prominent in H'^CO^ (and N2H^ not shown here) but not in 
C^'^S (and H2CO also not shown here) northwest of DR21(OH). 
In '^CO and '^CO 2— >1, this feature is discernible as well, 
but largely blended with other emission components. This sub- 
filament is the eastern end of the large-scale filament F3 that was 
observed in lower-angular resolution '^CO 1^0 channel maps 
(Fig. [3]l. It was also seen in isotopomeric CO 3— >2 data (Vallee 
& Fiege 2006) and interpreted by the authors as being caused 
by outflow emission from DR210H(M). An outflow is indeed 
detected at this position (see Sect. 4.3 1 but the emission in the 
velocity range — 1 to 1 km s"' is clearly due to an individual 
cloud fragment. Spectra towards a position within this fragment 
are shown in Fig. [8 An individual Gaussian line is seen at 
km s ' for all optically thin lines (at this position rather close to 
DR210H(M) also for C^'^S). It becomes obvious that the '-CO 
2^1 line - even at an angular resolution of 11" - is not useful 
for separating individual cornponents since it is too sensitive to 
low-density emission. Figure^superimposes NaH^ and H'^CO^ 
emission in the two major velocity ranges, the bulk emission 
between -6 to -3 km s"' in grey scale and the fragment emis- 
sion between -1 and 1 km s"' as overlaying contours. The emis- 
sion in H'^CO^ for both velocity ranges peaks at the position of 
DR210H-(M). A secondary weaker peak is found ~50" further 
soutii, close to N48 (or DR210H-(S)f] where N2H+ between -6 
and -3 km s ' has its peak emission. This supports the 'hot core' 
scenario for DR21-OH(M) because the gas there is warmer and 
chemically richer than in the dense, cold core of DR210H-(S). 
This core, however, is embedded in a warmer envelope traced by 
H'^CO^. Even more interesting is the km s ' component. The 
emission distribution in H'-'CO^ and N2H^ suggests that this 
fragment 'falls' on the densest part of the DR210H(M) clump. 
The N2H^ peak is more distinct and smaller in extent and shifted 
with respect to H'^CO^, again indicating that H'^CO^ emission 
arises from the warm envelope. 



4.3. Ou flows 

Figure [To] shows in more detail the outflow features detected in 
Fig.|6] From the position-velocity plot and positionally averaged 
'^CO spectra, we can clearly define the blue emission range be- 
tween -30 and -17 km s ' . The redshifted emission is more diffi- 
cult to discern because at higher positive velocities, a part of the 
line emission in the southern part of the DR21 filament is due to 
the 'Great Cygnus rift' (see Sect. 4.2.2 1. However, the velocity 
range ~17 to ~30 km s"' characterizes well the red wing. 

Several blue outflow sources are clearly recognized. These 
are all associated with mm-continuum sources, i.e. N53, N44, 
and N45. Source N51 shows broad line wings in '^CO and 



HCO^ (Fig. 14 1, which is most likely only the south extention of 
the N53 outflow. Southeast of N44 and at off'sets -100",-140" 
are clearly defined blue sources that have no red counterparts. 
Sources N53, N48, and N44 show outflow emission in SiO 2— >1 
(Motte et al. 120071) with a prominent red wing for N53 and a blue 
wing for N44. They are classified as massive infrared-quiet pro- 
tostellar cores (Motte et al. 2007) with masses of 85, 197, and 
446 M0, respectively. N45 exhibits no SiO emission but is asso- 
ciated with an H2 jet (source A 3-1 in Davis et al. 2007 ), which is 
probably driven by a low- or intermediate-mass YSO (all H2-jets 
are marked with crosses in the plot). N44 is the known outflow of 



* In this paper, we follow the notation of mm-continuum clumps 
given m Motte et al.l2007l This source list is more complete than the 
one from Chandler et al. 1993 However, for clarity, we still give the 
former names DR210H(Main, South, West) on occasion. 
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Fig. 9. Maps of fine integrated N2H+ 1^0 (top) and H'^CO+ 
1— >0 (bottom) emission of the DR21(OH) region. For N2H^, we 
integrate over all hyperfine structure components, the greyscale 
covering the velocity range -6 to -3 km s"' going from 0.8 to 
6.7 K km s ' and the black contours covering the range -1 to 
1 km s"' and go from 0.32 to 1.93 in steps of 0.32 K km s"'. 
For H'^^CO^ the same velocity range -6 to -3 km s ' represents 
grey scale (0.64 to 3.8 K km s"') and the black contours the 
velocity range -1 to 1 km s"' with contours from 0.67 to 6.6 by 
0.54 K km s"'). Green triangles mark mm-continuum sources, 
following the notation of Motte et al. ( 120071 ) . 



DR21(OH), seen in CS 5^4 emission (Richardson et al. [T994] ) 
and indirectly traced by the spatial and kinematic distribution 
of maser sources (Plambeck & Menten ,1990) . However, no H2 
emission was detected (Davis et al. 120071 ). 

This source is well studied (e.g. Garden et al. |I99IallI99Ibl 
' 1 9921 ). so we do not go into details here. In our map, we see also 
the classical east-west orientation of the flow with some north- 
south features (e.g. component 'C in Garden et al.). Among the 
other sources, only N53 and N44 show red wing emission. It is 
not clear what the emission feature running east-west at offsets 
-60" in declination is. In Fig. 10 it shows up in HCO^ emission 



as a typical outflow feature. But it has no blue counterpart and 
may well be associated with the Cygnus Rift. 
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Fig. 10. Map of line integrated ^^CO 2^ 1 emission in the DR21 
filament. Blue contours indicate a velocity range of -30 to -17 
km s-i (levels 12 to 60 by 8 K km s ' and then 100, 140 K km 
s"'), red contours a range of 17 to 30 km s ' (levels 6 to 36 by 
6 K km s"'). Green triangles indicate mm-continuum sources 
from Motte et al. (2007), black crosses mark H2 jets and knots 
discovered by Davis et al. (2007), and the red square south of 
DR21-OH indicates the 'Extremely Red Object' EROl detected 
by Marston et al. 620041) . 



5. Analysis 

This section is devoted to an analysis of the physical proper- 
ties of the DR21 filament. We, however, exclude the area around 
DR21 since this region is perturbed by outflow emission and line 
confusion. We focus instead on the northern part of the DR21 fil- 
ament, including the region of DR21(OH). 

5.1. Physical properties of tine sub-fiiaments 

Our molecular line observations show that the DR21 filament 
is connected to several sub-filaments seen on large and small 
scales. In the '^CO 1— >0 channel maps (Fig.|3]l, it becomes obvi- 
ous that the most prominent filament F3 covers a velocity range 
of ~3 km s"' (between -4 km s"' and -1 km s"') and has a pro- 
jected length of ~1 pc (the width is ~1.5 pc). The channel maps 
show that the top of the column lies at higher velocities, indi- 
cating that this part of the filament is tilted away from the ob- 
server. Assuming an average angle to the line-of-sight of 57.3° 
(in the case of a random distribution of orientation angles, the 
average angle is equal to this value), this implies a dynamic life- 
time of ~ 1.5x10^ yr The mass of the filament (for more de- 
tails see Appendix A about the mass determination from '^CO 
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Fig. 11. Greyscale map of line integrated NaH^ 1 ^0 emission 
of the DR21 filament with spectra of HCO^ overlaid in red. 
The spectra cover a velocity range of -10 to 10 km s ' and a 
temperature range of -1 to 12 K. Green triangles indicate mm- 
continuum sources from Motte et al. (2007) and are named ac- 
cordingly. N44, N38, and N48 coincide with the mm-sources 
DR210H-M(Main),W(West), and S(South), respectively, named 
by Chandler et al. (1993J. The 0",0" position is DR210H(M) 
(see Sect.[3]l. Spectra of the sources in bold (N53, N51, and N48) 
are displayed in Fig. 14 



data) is 2600 Mg and the average density is 690 cm"^. This sub- 
filament indeed corresponds to 'clump 7' in the DR21 region, 
identified in the '^CO 2^1 survey by Schneider et al. (120061) . 

The higher angular resolution IRAM observations resolve 
the detailed structure and show how this sub-filament connects 
to the DR21(OH) clump within the filament (see Sect. |4.2.3| l. 
The mass, determined from the line integrated H'^CO^ inten- 
sity (see Appendix A), is 1320 M© with an average density of 
8.1x10^* cm"-'. This mass estimate, however, is rather uncertain 
because of the uncertain [HCO^]/[H2] abundance ratio. From the 
channel maps (Fig. |7]i, we see that the filament extends up to ~ 
0.5 km s"' and therefore has a total radial velocity range of 4.5 
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Fig. 12. The grey area, divided by the total line-integrated inten- 
sity over all velocities (the 'normalized line intensity'), is shown 
in Fig. 13 Note that only spectra with a real infall line profile 



(T(blue)>T(red)) were considered. 
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Fig. 13. Greyscale map of the ratio of self-absorption to line- 
integrated HCO^ 1— >0 emission in the DR21 filament. Overlaid 
on that are red contours of H'-^CO^ 1— »0 emission (contours go 
from 3 to 13 Kkm s"' by 1.5 Kkm s"'). Green triangles indicate 
mm-continuum sources. The rectangle indicates the zoom region 
displayed in Fig.|9] 



km s"^ at a length of ~1 pc. If we assume that it connects to the 
large filament F3, the lower velocity levels are equally around 
-A km s"', thus implying a total velocity range of 4.5 km s"'. 
Again, assuming an average angle to the line-of-sight of 57.3°, 
this implies a very short dynamic lifetime of ~ 1.4x10^ yr. 

In addition, we note that in the high spatial resolution IRAM 
data, this sub-filament shows a position angle of ~45° in pro- 
jection to the sky, while on a larger scale F3 is mostly east-west 
oriented (PA ~ 100°). This bend, which could be related to the 
deep gravitational well of the DR21(OH) clump, may also be 
present along the line-of-sight, explaining the apparent accel- 
eration from 3 to 4.5 km s"' along the sub-filament between 
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Fig. 14. Spectra of optically thick ('^CO 2^1 and HCO+ 1^0) 
and thin (N2H+ l-»0, C^^-^S 2^1,and Hi3CO+) spectra from se- 
lected positions in the DR21 filament. The '^CO intensity was 
reduced by a factor of 4. To avoid confusion, the lower frequency 
HFS components of the N2H^ line were suppressed. 



large and small scales. Assuming a constant velocity of 4.5 km 
s"^ (the largest observed range of velocity) along the whole F3 
sub-filament, the angle of the sub-filament to the line-of-sight 
would be 48° before the bend and 90° after the bend, close to 
DR21(OH) (i.e. streaming away from the observer along the 
line-of-sight). This 3D view of F3 seems to be corroborated by 
the largest velocities in the IRAM data being found closest to 
DR21(OH). The total length of the F3 sub-filament would then 
be at least 9.4 pc (7 pc divided by sin(48°)) with no information 
about the length along the line-of-sight after the bend, leading 
to a total dynamic timescale of ~2.1xl0*' yr, and an input mass 
rate of ~1.9xl0""' Mq yr"' (total mass of F3 on large and small 
scales divided by the timescale). 



5.2. Infall signatures across the filament 
Figure 1 1 shows an overlay of HCO^ 1 - 



spectra (red) on a 
greyscale plot of line integrated NaH^ 1— >0 emission (both at 
an angular resolution of ~30")- It is obvious that HCO^ dis- 
plays along the whole filament the 'typical' profile of an op- 
tically thick spectral line which is indicative of infalling gas, 
i.e. a blue-shifted line wing that is stronger than the red wing 
(Myers et al. 1996). This self-absorption was akeady visible in 
'^CO emission (Fig.[6]l but becomes more pronounced using the 
HCO^ line. In order to quantify this infall signature, we deter- 
mined the normalized integrated intensity (see Fig. 12 between 



the blue and red line components (with T(blue)>T(red)), i.e. the 
'amount' of self-absorption, and plotted this quantity in Fig. [13] 
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It is remarkable that self-absorption is seen everywhere in the 
filament on a size scale of ~10'=5 pc with smooth transitions 
and not only locally at the position of mm-continuum sources. 
However, spectra with strongest self-absorption are found at the 
position of N44 (DR210H(M)). 

To exclude the possibility that we observe several line com- 
ponents, we display in Fig.[T4]examples of various lines at three 
selected positions (that are modeled in Sect. 5.3 1. First, the spec- 
tra show that '^CO is not the most reliable infall tracer since it 
is ubiquitous and the infall profile is most clearly seen in HCO^. 
The figure also clearly proves that not two different cloud com- 
ponents are observed since all optically thin lines (C^'^S, N2H^, 
H'^CO^) peak in the emission gap of HCO^. Apart from the 
known line components at -3 km s ' (bulk emision of the DR21 
filament) and +9 km s^' (W75N, only visible in '^CO 2-^1), 
the km s"' component (see Sect. 4.2.3 i is best traced in N2H^ 
emission. 

Since it is assured that we observe a true infall signature 
across the DR21 filament, we determined the infall velocity v,„ 
to first order using the method described in Myers et al. ( 1 19961 1 . 



1 + eiTBolTp) 



In 



(1) 



Here, cr is the H CO^ velocity dispersion, Trb, Vred and 
Tbd, ^hiue the temperature and the velocity of the red and 
blue part of the absorption profile, respectively, e=2.71828 
the Eulerian number, and To is the brightness temperature of 
the dip. In all calculations, we excluded the area containing 
DR21 {Ad <-100")- The velocity dispersion is determined 
from the FWHM line width of the optically thin H'^CO+ line 
(cr=FWHM(Hi3C0+)/V8 ln(2)). We disentangle the optically 
thin emission of the molecular ridge ( — 3 km s"') from the in- 
falling sub-filament at km s"' with a double Gaussian fit that 
precisely determines cr of the bulk emission to be -3 km s"' . For 
the optically thick HCO^ line, we did not perform such a double 
Gaussian fit since the -3 km s"' component, showing the infall 
signature, dominates strongly over the single Gaussian km s"' 
component (this means that the subfilament does not show an in- 
fall signature). We determined the temperature and the velocity 
of the red {Tkq, Vred) and blue (Tbd, Vbiue) part of the absorption 
profile by fitting a Gaussian to each wing. We only computed 
the infall speed in a region where the emission of the second 
component does not contribute significantly and confuse the ab- 
sorption profile. We note that we consider each point in the map 
as a single two-layer model (as outlined in Myers et al.11996'). 



Figure 15 shows the results of this procedure. Supersonic in- 
fall velocities, v,„ >0.2 km s"' for gas of a temperature of 10 K 
(typical value found), are detected within the red contour line, 
i.e. covering a large part of the filament. Interestingly, the high- 
est values of v,„ (~0.6-0.8 km s"') are not detected exactly to- 
wards the column density maxima (indicated by N2H^ and mm- 
continuum emission in the right panels of Fig. 15 1, but, instead, 
show some significant spatial offsets ranging from 20" to 30" 
(i.e. of the order of 0.2 pc) in projection. In the center region 
(DR210H(M) and N40) and in the south of the filament (N38 
and N48), high values of v,„ are mainly due to the increasing 
optical depth t of the HCO^ line since the infall velocity is 
determined by the ratio Tbd/Trd (eq. [TJ, which increases with 
T. In the north, higher values of v,„ are offset by ~20" to the 
dense clump containing N51 to N54 and are mainly related to 
the larger linewidths of the N2H^ line (see Fig. 18 showing the 
same increase in linewidth for H'^CO^ displaced with respect to 
the dense clump). 



Please note that this behaviour does not involve a gradual 
increase in infall speed since we only determine the infall veloc- 
ity at each point of the map assuming local infall at each point. 
From this figure, we can not extract the radial profile of the in- 
fall speed within the whole filament. For that, we need a more 
realistic radiative transfer model, by considering the density and 
temperature structure in the dense clumps of the filament as we 
present in the next section. 



5.3. Radiative transport modelling with Simline 

We used the ID non-LTE radiative transfer code Simline 
(Ossenkopf et al. 120011 ) to simultaneously model the observed 
optically thin H'^CO^ and optically thick HCO^ lines of the two 
clumps indicated on the right panel of Fig. 



15 These were ex- 



tracted from the mm-continuum survey obtained by Motte et al. 
(I2007I I and characterize well the two major regions of infall. The 
southern Clump- 14 is associated with DR21(OH) and contains 
the cores N36, N38, N41, N44, and N48. In the northern part of 
the filament, Clump-16 contains the cores N37, N43, N51, N53, 
and N54. The advantage of using mm-continuum data is that we 
have a rather reliable mass estimate than can be used as input to 
the model. Our ultimate goal is to derive the infall velocity and 
the level of depletion by constraining a simple physical model 
for these clumps. 



5.3.1. Input parameters and assumptions 

For both clumps, we used the Gaussian parameters (flux and 
FWHM) determined by Motte et al. (2007) and assumed the 
more realistic p(r) oc r^^ density profile to estimate the 90% 
mass and corresponding radius. We also recalculated the mass 
down to the level of 10% of the peak contour Clump-16 tak- 
ing a mass-averaged temperature of 20 K (rather than the me- 
dian temperature of Cygnus X clumps, 15 K, used by Motte et 
al. 2007), which is in closer agreement with NH3 measurements 
of this clump (Wienen 2008). The model self-consistently takes 
the mass corresponding to the mass-averaged temperature with 

a profile T(r) = r,„^^ j , where /3 is between -0.2 and -0.6 (T,,, 

and r,„ are the temperature and radius of the inner layer, see be- 
low). We adapted a two-layer geometry for the model since runs 
with only one layer did not reproduce well the observed line pro- 
files, a result similar to that of Jakob et al. (2007) who equally 
applied a two-layer model for their low-angular resolution mid-J 
CO and atomic CI lines. However, they modeled their observa- 
tions without considering molecular depletion and no or only 
very low infall speeds. We applied a subthermal external layer at 
constant temperature with an average density of n^lO-'-lO'* cm"-* 
because the self-absorption dip is very deep. 

For the determination of the turbulent line-width, we used 
the optically thin line of H^^CO^. A single Gaussian fit yielded 
line-widths of ~2.5±0.5 km s"' for Clump-14 and ~1.5±0.5 
km s ' for Clump-16, respectively. This line-width is a result 
of thermal width, internal turbulence, and possible infall. Since 
both turbulence and infall broaden the line, their parameter space 
must be explored carefully. 

For the HCO+ abundance, we ta ke a value of X(HCO^) 1- 
5x10"' (e.g. Marseille et al. 120081) and adopt an isotope ra- 
tio i^c/'^C = 67 (e.g. Lucas & Liszt 1998). However, HCO+ 
is known to be depleted in cold environments and at densities 
higher than ~10^ cm"^ (Tafalla et al. 2002 ). We expect this hap- 
pens in Clump-16 but to a lesser extent in Clump-14 since this 
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Fig. 15. Left: Color plot of infall velocity. The red contour line indicates the thermal sound speed of 0.2 km s White triangles 



mark mm-continuum sources. The two red lines indicate the cuts where we model spectra (Sect. 5.3 i. The cuts were selected to cross 
the center of each ellipse. Middle: Contours of infall velocity (levels are 0.3, 0.5, 0.7 km s"') overlaid on a plot of line-integrated 
N2H^ intensity determined from a Gaussian line fit to the -3 km s"' component. It becomes obvious that the highest infall velocities 
are offset from density maxima. Right: Mm-continuum emission map taken from Motte et al (120071 ). where the two ellipses indicate 
the clumps considered for modeling and calculating infall properties. Note that all plots cover exactly the same area and have the 
same scales. 



region contains a hot core. We thus introduce a threshold density 
niresh in our model, with a sharp jump in the abundance profile 
if n> n,„i/,. The depletion factor is assumed to be in the range 
1-100. 

Our strategy to constrain the infall speed was to vary the de- 
pletion ratio, the abundance, and the threshold density (between 
the given limits). To determine the best fit to the model, we per- 
formed a reduced x^-is&l- 



5.3.2. Results 



Figures 16 and 17 show the observed and modeled spectra of 
HCO+ and H'^CO+ for Clump-14 and Clump-ie, respectively. 
For better orientation, the two cuts are indicated in Fig. [TS] 
Table |2] summarizes the explored parameter space and gives the 
results for the best fitting model. We note that the spectra of 
both sources have broad wings caused by outflow emission that 
we did not account for in our modeling. 





Clump-14 


Clump- 16 


total mass 


4900 Mo 


3346 Mo 


size 


0.52 pc 


0.615 pc 


density profile (0-2) 


-8.0 


-8.0 


mass-averaged T 


20 K 


15 K 




-0.6 


-0.5 


line width 


2.5 km S-' 


1.5kms-' 


X(HCO+) 


2.5x10-' 


2.5x10-' 


'2C/"C 


67 


67 


depletion ratio 


7 


28 


Thco+ 


25 


26 


'''h"co+ 


0.81 


0.99 




-0.6 km s ' 


-0.5 km s- 



Table 2. Results of Simline modelling where 0-2 is the density 
profile exponent at the transition to the subthermal layer and 
X(HCO+ the HCO+ abundance. 



Clump-14 

For Clump-14, our best-fit model closely reproduces the spectral 
lines of HCO+ (r ~ 25) and H'^CO+ (t ~ 0.8). It was necessaiy 
to use a rather large outer layer to reproduce well the absorption 
dip, which extends to R ~ 0.9 pc increasing by 20% the mass of 

our model. We found that a temperature profile of p 0.6 and 

a mass-averaged temperature of 20 K provides a closest fit to our 
data, where the discrepancy in the peak intensity of the lines is 
smaller than 20%. This also indicates that the inner part of this 
clump tends to be hotter, reaching up to ~100 K in the central 
part. 

The depletion occurs at a radius of 0.25 pc with a depletion 
ratio of ~ 7. This may be due to the temperature because the 
central region is a hot molecular core. With this configuration 
of geometry and parameters, we derive an infall speed of — 0.6 
km s ', which compares well with the values derived from the 
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Fig. 16. Observed (black) and modeled (red) HCO^ and H'-'CO^ 1— >0 spectra of Clump-14 (see Fig 15 for the location of the cut). 
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Fig. 17. The same line Fig. 16 but for Clump- 16. 





Myers et al. method (see Sect. |5.2[ ). 
Clump-16 

The best-fit model for Clump-16 (Fig. [TTjl reproduces very well 
the H^^CO+ line (t ~ 1) but less well the optically thick HCO+ 
line (t ~ 26). Since the densities are lower for this clump, the ex- 
tended outer layer is smaller than that of Clump-14 and reaches 
out to Rj i, ~0.8 pc reproducing the 10^^ cm"^ density. As we kept 
the exponent /3 for the temperature profile free, we found that 
j3 = -0.5 most closely fits our data. This results in a maximum 
of T~60 K in the central region. The depletion occurs at R~ 0.2 
pc with a depletion ratio of 28, implying that in the central, very 
dense part of this clump, this molecule is strongly depleted. 

The fit value of the infall speed is — 0.5 km s ' , again within 
the range obtained from the more simple Myers et al. method. 
The strongest constraint on the infallvelocity is the ratio of the 
left-to-right wing intensity ratio, because the position of the plot- 
ted spectra is aligned with the optically thin line. 



Table 3. Infall properties for Clump 14 and 16, i.e. the two main infall 
regions. 





R" 




<nH, >^ 


V'' 

111 


dM/df 




[pc] 


[Mo] 


10^ [cm-3] 


[kms-'] 


[Moyr-'] 


Clump 14 


0.52 


4900 


1.26 


0.6 


5.8 10-" 


Clump 16 


0.615 


3346 


0.52 


0.5 


2.8 10-3 



" Equivalent radius of infall region 

* Mass from mm-continuum (see Table 2) 
average H2 density, derived from the mass given in col umn 2 
infall speed, determined from Simline modelling (Sec. 



mass infall rate M « M/t = M v,„/R 



5.3 1 



5.4. Global infall properties of the DR21 filament 

Table 3 summarizes the infall properties of the two main infall 
regions that we defined in the Simline modeling. It is first im- 
portant to understand that we discuss here the global dynamic 
features of the molecular clumps on a size scale of ~0.5 pc, and 
not the physical properties of individual cores on a size scale of 
-0.1 pc. 

The large global infall speeds on a ~0.5 pc size scale and 
the high masses of the clumps lead to large global infall rates 
of M of 5.8x10-^ Mo yr-i and 2.8x10"^ Mq yr"' for Clump- 
14 and Clump-16, respectively. However, one has to distinguish 
between large-scale flow motions and local core collapse. The 
expression 'infall' , representing material that is falling onto a 
protostar or - more generally - into a gravitational potential 
well that could form many stars within a core that is gravitation- 
ally collapsing has to be used with care. In Fig. 15 we see that 
the highest values of 'infall' speed do not always correspond to 
dense cores, only source N48 shows a local peak for v,„. For 
the rest of the DR21 filament, the highest 'infall' speeds are 
found in the less dense western part of the filament and not in 
the dense post-shock gas. In view of the high mass infall rates, 
this implies that there is a significant contribution of matter in- 
put either by directed flows (as would be the case for large-scale 
convergent flows) or gravitational collapse of the DR21 filament 
itself on a large scale. Both clumps have a low virial parame- 
ter a[^of 0.46 and 0.30 for Clump 14 and Clump 16, respec- 
tively, and are thus strongly gravitationally bound. Their free- 
fall Hfetimes t/y|^are 0.9x10^ yr {Clump 14) and 1.3x10^ yr 
{Clump 16), a factor 20-30 shorter than the sound crossing time, 
which indicates that they are indeed in global collapse (if the 
simple picture of a collapse onto a single gravitational potential 



' The virial parameter is defined to be a=Myir/M where the virial 
mass is M,,,,[Mo]=699 cr^[km/s]"^ R [pc] and the mass M from Table 2. 



tyy[yr]=1.37 10' 



fi( Whrrrl) 
\ 2,KH2) / 
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is maintained). On the other hand, the infall rates are in agree- 
ment with the theoretical predictions of both the micro-turbulent 
core model (McKee & Tan 2003 ) and the gravo-turbulent frag- 
mentation model (e.g. Klessen et al. 1 19981 [20001120011 Padoan 
& Nordlund l2002l Banerjee & Pudritz l2007l l. 

5.5. The kinematic structure of the DR21 filament 

To further characterize the kinematic structure of the DR21 fil- 
ament, we determine the line velocity and velocity dispersion 
using H'^CO^ and NaH^ by fitting a Gauss ian profile to the -3 
km s ' component as described in Sect. 5.2 For N2H^, we addi- 



tionally fitted the hyperfine structure pattern of the molecule. 



5.5.1. Velocity dispersion 



Figure 1 8 reveals that the H'-'CO^ main beam brightness temper- 
ature map shows a V-shaped feature in the north of the filament. 
The western component is a factor 2-3 weaker than the eastern 
one, which has a pronounced peak at the position of N5 1 . This 
peak is stronger than that towards N44/DR210H(M) (the reverse 
is true for the hne integrated intensity, see Fig. [5]). 



H cn 1 D line temperature 



o - FWHM/(B ln(2))°- 





50 ' too 



Fig. 18. Color scale maps of the results of Gaussian line fitting 
(temperature and velocity dispersion) to the -3 km s"' compo- 
nent of H'^CO+ 1^0 along the DR21 filament (but excluding 
the region of DR21). Contours of velocity dispersion (second 
panel) at 0.75 and 1 km s ' are overlaid on both maps. 

Most remarkable, however, is the velocity disper- 
sion cr (determined from the FWHM of the line with 
cr=FWHM/(81n(2))°^) seen in H^^CO^ as well as in N2H+ 
(not shown here). The velocity dispersion increases towards 
the geometrical center of the filament from sonic levels to 
values of ~1.3 km s In the northern part of the filament 
(declination offset>0"), the gradient is stronger than in the 
southern part (south of N44/DR210H-(M)) where we see 
a more homogeneous distribution of cr. Interestingly, in the 
northern part of the region highest cr-values are found between 
the two vertical features of H'-^CO^ emission (with a decrease 
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360 Fig. 19. Color scale maps of the results of Gaussian line fitting 
(line integrated intensity and line velocity) to the -3 km s"' com- 
ponent of N2H^ 1— >0 along the DR21 filament (but excluding 

^'^ the region of DR2 1 ) . Contours of velocity dispersion of H ' ^ CO^ 
(0.75 and 1 km s"') are overlaid on the area map. Contour lines 

jj,) of the line integrated (-1 to 1 km s"') emission of the km s"' 
component are overlaid on the velocity map (levels 1, 2, 3 K km 



of cr from ~1 km s ' to ~0.8 km s ' in the densest region in 
the east). This becomes even more obvious in Fig. 19 where 

N2H+ 



we plot contours of cr (from H'^CO^) over a map of 
line integrated intensity (proportional to the column density). 
Thus, there is a clear offset between high cr and high column 
densities. This scenario can be explained if it is assumed that 
clumps are produced by turbulent flows because in this case, 
the densest gas has been shocked and slowed down so that the 
largest velocity dispersions do not occur in the densest regions 
but in the outskirts (Klessen et al. 12005 . Vazquez-Semadeni 
et al. 12^8 Feden-ath et al. I2009I) . We return to this point in 
the discussion. The southern part (offset DEC<0") seems to 
behave differently. The N2H^ column density and H'^CO^ line 
temperature are highest in the core region (~0.5 pc diameter) of 
N38, while the velocity dispersion remains on a rather constant, 
high-level of around 1.2 km s ' in an area of ~1 pc diameter. 
However, a projection effect could explain the difference.: it is 
possible that a similar offset between velocity dispersion and 
highest density region is seen face-on. 

5.5.2. Velocity field 

Figure[T9]shows the complex velocity field of the DR21 filament 
determined from N2H^ (the H'^CO^ map not shown here looks 
very similar). Three main, basically horizontal, projected veloc- 
ity gradients are observed. While the northern and southern part 
show a decrease of velocity from west to east (-2.2 and -0.8 km 
s"' pc"'), the middle part has a positive value of -1-2.3 km s"' 
pc"'. Interestingly, we observe in this region a 'turnover' in the 
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velocity gradient that corresponds to the material of the km s"' 
component 'falling' onto the bulk of the filament (see Fig. 19 1. It 
is clear that the velocity pattern cannot be explained with a sin- 
gle rotation of the filament along a vertical axis. An alternative 
explanation would again be that at least a part of the observed 
motion is due to convergent flows (see Sect.|6]l and/or rotation of 
the individual parts of the filament. 

6. Discussion 

The dynamic structure of the DR2 1 filament exhibits several re- 
markable features: 

- Several sub-filaments are connected to the DR21 filament, 
the most massive one, F3, being clearly linked to the 
DR21(0H) clump. 

- Global infall signatures are seen all across the filament (on a 
length scale of a few pc). 

- The highest values of the velocity dispersion are found in a 
vertical projected region of low column density gas, offset 
from the dense clumps seen in N2H^. 

- Three horizontal (west-east) velocity gradients (-0.8, -h2.2, 
+2.3 km s"' pc"') exist in the filament. 

All of these findings point towards a very dynamic nature 
of the DR21 filament in which motions are not randomly dis- 
tributed but show clear signs of systematic dynamics. These fea- 
tures all point to local converging flows being the driving source 
of the whole filament, with an increase in density at the stagna- 
tion points of the flows (Sect. |6.1| l. In addition, the large-scale 
infall signatures observed all over the filament point towards a 
scenario in which gravity is the dominant driving force on the 
scale of the filament (Sects. [6^ and |6.3| l. On the other hand, the 
origin of the large-scale flows is difficult to constrain. It is possi- 
ble that the large-scale dynamics are a remains of the formation 
process of the molecular complex, but may also originate either 
from self-gravity on larger scales, or large-scale turbulence in- 
jection (Sects. 6.4 and 6.5 1. 



6.1. Organized motions: signposts of convergent flows 

The clouds of the whole Cygnus North region are very similar 
in their physical properties (morphology, density etc.) to clouds 
that were created by colliding flow scenarios in models (e.g. the 
final stage of the Gf2 model from Heitsch & Hartmann 2008). 
The DR21 filament is part of a complex network of intercon- 
nected filamentary structures (Sec. 



4.1 



and could thus indeed 
be the result of a convergent, turbulent flow. In particular, the 
offset of the infall speed from the densest regions is consistent 
with a scenario where supersonic flows shock at their stagnation 
points are offset from the dense postshock gas on both sides of 
the filament. In this case, high-density clumps/cores are built up 
at the stagnation points of the colliding flow and would explain 
the distribution of the dense mm-cores along the DR21 filament. 
It suggests that the cores are dynamically distributed through the 
region before they evolve into intermediate or high-mass star- 
forming clusters. From the mm-continuum survey of Motte et 
al. (2007), we determined that approximately 15% of the mass 
of the filament is contained in dense cores, in accordance with 
the predictions of Heitsch & Hartmann (I2008I I. 

Another observational results that agrees with supersonic 
isothermal simulations is the velocity field and the distribution 
of the velocity dispersion (cr). We observe different velocity gra- 
dients along the filament. Figure [19] shows in the northern and 



southern part of the filament a negative gradient in the WE- 
direction whereas at the position of DR21(OH) (where the sub- 
filament falls onto the densest region), this gradient is reversed. 
It is unlikely that the gradients are only caused by rotation (see, 
e.g. Goodman et al. |1993| for typical observational signatures of 
uniform rotation). Moreover, they can be explained in the tur- 
bulent colliding flow view, where clumps are sheared and com- 
pressed due to the flow motions and retain the signature of the 
external flow that formed them (Ballesteros-Paredes et al. 119991 
Vazquez-Semadeni et al. i2008i) so that only parts of the filaments 
are rotating. 

The shock front created by the colliding flows creates an in- 
ner dense post-shock region and an outer lower-density region 
of supersonic inflow. Cores are then formed by 'gravoturbu- 
lent fragmentation' (Klessen et al. '2005' Hennebelle & Chabrier 
2008), which characterizes a hierarchy of smaller and smaller 
entities. The local collapse of a core occurs when self-gravity 
overcomes the gas pressure (including the turbulent contribu- 
tion to pressure). The whole process leads to a higher velocity 
dispersion outside the clump/core than in the more quiescent 
post-shocked interior (Klessen et al. 2005 Gomez et al. 20071 
Vazquez-Semadeni et al. I2008I I. which is precisely what we ob- 
serve in the DR21 filament. Figure 1 in Klessen et al. ( 12005 l l il- 
lustrates this scenario in a non-magnetic SPH (smoothed particle 
hydrodynamics) model. Here, we provide evidence of the same 



scenario in Fig. 20 where we show a map of the velocity dis- 
persion obtained in a driven turbulence simulation by Federrath 
et al. (i2009i) . The simulation used fully compressive (curl-free) 
forcing to excite turbulent motions, though models with solen- 
odial (divergence-free) forcing (not shown here) produce similar 
results. In general, compressive forcing produces greater den- 
sity contrasts with higher density enhancements and larger voids 
than solenodial forcing (Federrath et al. |2008l|2009l l. This could 
be the case for the DR21 filament since it is here where we ob- 
serve the massive dense cores with densities higher than 10^ 
cm"^ and masses up to 200 Mq. The velocity field primarily con- 
sists of converging flows, which create strong density enhance- 
ments. Column density contours are overlaid on the velocity dis- 
persion, showing that cr and density peaks are offset from one 
another, as in our observations shown in Fig. 15 (middle panel) 
of the DR21 filament. High density gas is primarily found in 
regions of low cr, while lower density regions have higher cr. 
We emphasise that the simulation did not include self-gravity. 
Nevertheless, the velocity dispersion-density structure is similar 
to what is observed in the DR21 filament, which indicates that 
the observed structures can be produced by converging super- 
sonic gas flows only, without the necessity of a global gravita- 
tional collapse. However, converging gas flows are also created 
by global collapse. If supersonic flows and global collapse act 
simultaneously, it is difficult to disentangle their individual con- 
tribution to the creation of the velocity dispersion-density struc- 
ture seen in our observations and in the simulations (see Sect. |6.5| 
below for a general discussion of the origin of converging flows). 

6.2. Self-gravity, the main driving force in the DR21 filament? 

Observationally, the ubiquity of the blue-shifted asymmetric line 
profiles in all optically thick lines indicates inflow motions. This 
is a strong hint that self-gravitation plays a significant role for the 
dense and massive DR21 filament. In addition, the most massive 
sub-filament F3, which is a good candidate to flow towards the 
main filament, seems to bend towards the center of the largest 
gravitation well of the whole filament, the DR21(OH) clump 
(see Sect. |5.1| l. This bend may be indicative of a dominant role of 
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Fig. 20. Results of a hydrodynamic turbulence simulation by Federrath et al. ( 120091 ). considering only compressive (curl-free) forc- 
ing. The Mach numbers reached in this simulation are comparable to the typical Mach numbers observed in the DR21 filament. The 
color-coded velocity dispersion is in units of the sound speed, while contours of column density are overlaid on it. The simulation 



was performed with the grid code FLASH3 and used 1024^ computational elements. This figure can be directly compared to Fig. [15 
(middle panel), also showing a prominent offset between the maxima of velocities and column densities. 
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Fig. 21. In the left-hand plot, we show the column density distribution from an MHD code (Hennebelle et al., Tevssier 120021 
Fromang et al. l2006l l. In the right plot, the projected velocity field is overlaid on the column density distribution. 



gravity as the flowing material approaches the main filament. At 
larger distances, the effect of large-scale magnetic fields could 
be the cause that keeps it in an almost east-west direction (see 
Sect.lOl. 



As shown in Sect. 6.1 



we are able to explain the spatial dis- 
placement between the maximum of the velocity dispersion and 
the dense gas and the velocity field of the DR2 1 filament in the 
context of supersonic isothermal simulations that do not include 
self-gravity. However, models of isothermal, non-magnetic, but 
self-gravitating turbulence produce a velocity divergence of typ- 
ically 0.6 km s ' pc ' for a 1 pc region (Vazquez-Semadeni et al. 
I2008I I. The authors propose that the assumption that the veloc- 
ity dispersion is produced only by random motion is incomplete 



and that even in the presence of driven turbulence a part of the 
observed velocity dispersion is caused by clump-scale inward 
motions due to gravitational collapse. In this large-scale inflow 
(LSI) scenario, the velocity field on all scales includes a signifi- 
cant inflow component. 

6.3. Role of magnetic field: MHD with self -gravity modeling 

We observed several sub-filaments that are attached to the DR21 
filament (Sects. |4.2| l. The orientation of the sub-filaments is 
not arbitrary: they follow the magnetic field lines that run or- 
thogonal to the NS-orientated DR21 filament (Vallee & Fiege 
120061 Kirby 120091 . Magnetohydrodynamic turbulence models 
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Fig. 22. Results from the same MHD turbulence model (Hennebelle et al., Teyssier'2002' Fromang et al. I2006I I shown in Fig. 20. 
Top: Velocity dispersion (left) and velocity field (right). Bottom: Magnetic field with field vectors overlaid on it. 



precisely produce such filamentary structures (Fig. 21 1 that are 
aligned with the magnetic field (Fig. [22] l. Hennebelle & Audit 
( 12008b found that self-gravitating MHD models produce more 
filamentary structures than purely hydrodynamic models. These 
numerical simulations were specifically designed for this object 
and the initial conditions resemble those used in Peretto et al. 
(120071) . They consist of an elongated clump with an initial aspect 
ratio of 2 and a density profile p(r,z) - P()/(l +i''/fo)'^ + (z/zo)^), 
where r = -s/ix^ + y^), zo = 2ro, ro = 5 pc and po - 500 cm"-*. 
The density at the edge of the cloud is equal to 50 cm"^ and is 
10 times lower than the value outside the cloud. Turbulence is 
seeded initially in such a way that the clump is approximately in 
virial equilibrium. The initial temperature is 10 K and the clump 
is threaded by a magnetic field parallel to the x-axis. Its inten- 
sity is proportional to the cloud column density and the peak 
value is about 7 jiG making the magnetic energy about 4 times 
lower than the turbulent one and 5 times higher than the thermal 
energy. This simple initial configuration for the magnetic field 
is suggested by the observations of Vallee & Fiege (120061) . who 
show that the magnetic field is indeed perpendicular to the major 
axis of DR21. The simulation was performed with the Ramses 



code (Tevssier 120021 Fromang et al. 120061 1. Since the cloud oc- 
cupies a small fraction (about 5 percent) of the computational 
domain, the following strategy is adopted. Initially a uniform 
128^ computational grid is used (level 7). Then, all the cells of 
density above 10 cm"^ and 40 cm"^ are refined to the level 8 and 
9, respectively. This ensures that initially the clump is described 
with about 3x10'' computing cells. As the collapse proceeds, the 
resolution is adjusted requiring at least 10 cells per Jeans length 
up to level 14. 

In the left-hand panel of Fig. |2T| the distribution of the col- 
umn density is shown and in the right-hand figure panel, the 
velocity field vectors are overlaid on top. The column den- 
sity is smoothed to a beam of ~0. 1 pc angular resolution (that 
of the IRAM 30m molecular line observations). The resem- 
blance to the DR21 filament is evident with an elongated verti- 
cal main filament and several subfilaments attached to it. Several 
clumps/cores (size scales <0.5 pc) with high H2 column densi- 
ties (>10^"* cm"^) have formed in the filament, in agreement with 
our observations. 

Figure 22] shows the mean velocity dispersion, the mean ve- 
locity field along the line of sight, and the averaged magnetic 
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field. For the velocity dispersion cr, we observe the same offset 
between regions of high density and cr as seen in Fig. |20] The 
modeled velocity field is as inhomogenous as the observed one, 
though the gradients have slightly larger values (up to a few km 
s"' pc"') and the same direction (no turn in sign as we found for 
the region around the DR21(OH) clump). The magnetic field is 
particularly interesting. The total magnetic intensity is displayed 
in Fig. 



22 



; (bottom), the arrows showing the direction of the mag- 
netic field in the xz plane. The sub-filaments are clearly aligned 
along the magnetic field and even sometimes pinched by it. This 
is in full accordance with our observations and points toward a 
scenario in which the DR21 filament is fed with material along 
EW- orientated magnetic field lines. A full network of these sub- 
filaments was detected in recent JCMT '^CO 3— >2 data (Richer 
et al., in prep.). A more detailed discussion of these observations 
and MHD modelling will be presented in a forthcoming paper 
(Csengeri, Hennebelle et al., in prep.). 

However, the prinicipal driver of these flows remains unclear 
Is it external, large-scale driven turbulence or the gravitational 
attraction of the gas itself, or both? A possible scenario is that 
first, the turbulence creates the density fluctuations (clumps and 
filaments) by means of locally convergent flows, and then, for 
some of the more massive fluctuations,gravity takes over, mak- 
ing the inflow stronger and driving the material to higher densi- 
ties. We probably see both effects in the DR21 filament: turbu- 
lent compression and self-gravity. The northern part is probably 
in an earlier stage where the external flow is still important in 
shaping the velocity field (hence the offset of the 1st and 2nd 
velocity moments from the density map). In the southern part, 
we see a more advanced stage in which gravity dominates and 
gravitational infall extends all the way to the center (leading to 
the coincidence of the density and velocity structures). 



6.4. Are the observed sub-filaments able to maintain the 
formation of the DR21 filament? 

We observed that several sub-filaments are attached to the DR21 
possibly serving as a reservoir of gas for 



filament (Sect. 4.2 



further mass growth. These organized flows help to sustain high 
accretion rates for times long enough to build up very massive 
clumps/cores (Vazquez-Semadeni et al. |2009l l. We note that the 
DR21(OH)-clump has a mass of higher than 5000 Mq in a ra- 
dius of ~0.5 pc, which makes it one of the most massive clumps 
known in the Galaxy. The F3-filament (Sect. 5. 1 1 has an average 



mass input rate of 1.9x10" M© yr" so that it may contribute 
up to ~ 30% of the total mass infaU rate of 5.8x10"^ Mq yr"' 
observed for the DR21(OH) clump as a whole (Clump 14, see 
Table 3) for the next 2.1x10* yr. Accounting for the other less 
massive sub-filaments, which may also contribute, this indicates 
that the mass flow rates are of the correct order of magnitude to 
explain the formation of the DR21 filament. On the other hand, 
they seem insufficient to build the present 5000 Mq DR21(0H) 
clump, which suggests that the mass flow rate must have been 
higher in the past and that the process that produced the DR21 
filament may be declining. 



6.5. Origin of the converging flows that formed the DR21 
filament 

One physical scenario that intrinsically includes strong dynam- 
ics is that of molecular clouds and complexes being formed by 
large-scale HI convergent flows, which would then drive the 
whole dynamics of the complex down to the star-formation scale 



(Kles sen et al. [T998l Foli ni & W alder l2006l Va zquez- Semadeni 
et allmT. Heitsch et al. l2008l Banerjee et al. 12009' Federrath 
et al. 2009). Common for all 'convergent-flow' scenarios is the 
basic principle that the clouds result from a pileup of material, 
i.e. atomic hydrogen, from large-scale supersonic flows, caused 
by energetic events such as supernovae, expanding Hll-regions, 
galactic spiral density waves. The driving force could be the in- 
jection of supersonic turbulence on large scales in a pre-existing 
molecular complex. Supersonic turbulence necessarily leads to 
organized supersonic motions on smaller scales, and then leads 
to the formation of filaments and clumps, and then stars. Finally, 
the self-gravity of a pre-existing molecular complex may play a 
role in driving the flows. Altogether, it is difficult to trace back 
the origin of the large-scale motions observed here. But at least 
on the observed scales, the motions are supersonic and are con- 
verging to form dense structures such as the DR21 filament. 
Whether a more quasi-static evolution is possible on smaller 
scales or in other parts of the complex to drive the evolution 
towards (high-mass) star formation cannot be addressed by the 
present observations. In Csengeri et al. (120101 ). we address the 
kinematic state of the small scale, densest structures in the com- 
plex, the i.e., massive dense cores, using PdBI data (Csengeri et 
al. in prep). 

To distinguish the formation scenarios more clarly on large 
scales, it is not clear which other observational signatures can 
be expected. The DR21 filament itself is already a rather dense 
structure and in order to detect residues of the formation pro- 
cess, it is neccessary to inspect the distribution of low-density, 
subsonic gas using atomic hydrogen and low-density tracers 
such as '^CO. We will address this point in a subsequent paper 
(Csengeri, Hennebelle et al., in prep.). 

7. Summary and conclusions 

We have presented a detailed molecular line study of the 
molecular ridge containing the star-forming regions DR2 1 and 
DR21(OH). This ridge is embedded in a large-scale network of 
filamentary structures, revealed by our maps of '■'CO 1— >0, CS 
2^1, and N2H+ 1^0 emission obtained with the FCRAO. It 
is the most massive (around 30 000 Mq) and dense (average 
density ^10"* cm"^) filament within the region and is labeled 
by us the 'DR21 filament'. Several sub-filaments are linked to 
the DR21 filament, the most massive one runs orthogonal to the 
North-South oriented ridge and has a mass (determined from 
'^CO 1^0) of 2600 Mq and an average density of 690 cm ^ 
Its inferred dynamical time is ~2xl0* yr. 

Higher angular resolution IRAM molecular line observa- 
tions in HCO+, H'^CO^, i^CO/i'CO 2^1, C^^'S, N2H+, and 
H2CO resolve the detailed structure of the DR21 filalment. The 
H'^CO^ 1— >0 data show how the sub-filament seen in '■'CO con- 
nects direcdy to the DR21(OH) clump. From the '^CO 2^1 
line mapping, we confirmed the known outflow sources DR21 
and DR21(OH) and detected three new ones, correlated with the 
mm-continuum sources N53, N44, and N45 (Motte et al. 2007]). 
The HCO^ 1— >0 line shows self-absorbed lines across the whole 
filament. Since optically thin lines peak in the gap and the blue 
wing of HCO^ is more intense than the red one, we conclude 
that this emission feature is due to infalling gas. The typical in- 
fall speed, determined with a simple method described in Myers 
et al. (119961) . is 0.6-0.8 km s"'. A more sophisicated non-LTE 
modelling of the HCO^ and H'^CO^ lines using the Simline ra- 
diative transfer code yields an infall speed of 0.5 and 0.6 km 
s ' for the northern and southern part of the DR21 filament (but 
excluding the DR21 region itself), respectively. 
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The kinematic structure of the DR21 filament is remarkable. 
We measured, using the N2H+, H'^CO^, and HCO^ maps, the 
highest values of the velocity dispersion in a vertical column of 
low column density gas that is offset from the dense clumps seen 
in NaH^. These results can be explained if the filament was pro- 
duced by turbulent flows. In this case, the densest gas has been 
shocked and slowed down so that the largest velocity dispersions 
do occur not in the densest regions but in the outskirts (Klessen 
et al. 120051 Vazquez-Semadeni et al. I2008I I. We also observe 
three horizontal (west-east) velocity gradients (-0.8, -2.2, -1-2.3 
km s"' pc"^) in the filament in which the one with the positive 
value marks the location where the filament is 'falling' onto the 
DR21(OH) clump. This velocity pattern cannot be explained by 
a single rotation of the filament along a vertical axis. An alter- 
native explanation would be that at least part of the observed 
motions is due to convergent flows. 

By comparing our observations of the DR21 filament 
with a hydrodynamic (Federrath et al. 12009 l l and a magneto- 
hydrodynamic turbulence model (Heimebelle et al., Teyssier 



12002 , Fromang et al. 2006 ), we infer that a very dynamic and 
fast mode of star formation occurs in the filament, in which gas 
is continuously replenished by subfilaments attached to the main 
filament. These subfilaments are aligned with the magnetic field 
direction that is perpendicular to the DR21 filament (Vallee & 
Fiege |2006l l. The DR21 filament is globally collapsing. All our 
observational findings are incompatible with the view of a quasi- 
static, pressure-bounded clump scenario. 
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Appendix A: Calculation of physical parameters 

A.1. Abundances 



We use a standard [^^CO\l[^^CO\ ratio of 67 (Lucas et al. [T998] l. 

The [HCO^]/[H2] abundance in high-mass star-forming re- 
gions is typically quoted to be 1.5-2.5 x 10"^ (Peretto et 
al. L2006, Marseille et al. 2008 ). For our mass determination 
from H'^CO^, we apply the a bundance value obtained from the 
Simline modelling (Sect. 5.3 1, which is 2.5 x 10"^. 



A.2. Determination of mass from molecular lines 

The total column density of any optically thin molecule can be 
determined from 

N[cm-^] = /(r„) j T^blKjdvlkms-'] 



(A.1) 



where 



f(Te.) 



3hZ 



and the Partition function is given by 



Z = 



hv 



+ 1/3 



(A.2) 



(A.3) 



in which h and k denote the Planck and the Boltzman constants, 
respectively, v is the frequency [GHz], ju is the dipole moment 
[Debye], J, is the upper value of the rotational quantum number 
and j Tinbdv is the velocity integrated line intensity on a main 
beam temperature scale 



hv 



kiexp(hv/(kT,,) - 1) 
and J(Tbc) = J(2.7K). 



(A.4) 



The values for the H'^CO+ 1^0 transition are v=86.75433 
GHz, hv/k=4.163 K, yu=3.3 Debye, and J{2.7K)=1.13. 

The temperature-dependent factor f(Te^) is (2.43, 2.84, 4.39, 
9.4)xl0i2 for 10, 20, 40, and 100 K, respectively. 

The values for the '^CO 1— >0 transition are 
v=110.201 GHz, hv/k=5.29 K, ^'=0.112 Debye, and 
7(2.7/r)=0.868. The temperature-dependent factor f(T(,i) is 
(0.91,1.25,1.66,2.08)xl0i5 for 10, 20, 30, and 40 K, respec- 
tively. 

The H2 column density from H'^CO^ is calculated to be 



N(H2)[cm"^] = 2.43 x 10'°N(H'^CO+) 



and for '^CO 1^0 from 
N(H2)[cm"2] = 4.7 x lO^N('^CO) 



The final mass is determined to be 
M[Mo] = N(n2) d^A(7T/mf2m„fi, 



(A.5) 



(A.6) 



(A.7) 



87rVV, [1 - exp(-hv/kT,MJ(Ted - J(Tbc)) 



where the projected area A is in square degrees, the distance d in 
parsec, the mass of hydrogen m//, and fi-1.36. 
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